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Summary. — With the start-up of the LHC, we can hope to find evidences for 
new physics beyond the Standard Model, and particle candidates for dark matter. 
Determining the parameters of the full underlying theory will be a long process 
requiring the combination of LHC and ILC data, flavor physics constraints, and 
cosmological observations. However, the Very Early Universe properties, from which 
the relic particles originate, are poorly known, and the relic density calculation can 
be easily falsified by hidden processes. We consider supersymmetry and show that 
determining the underlying particle physics parameters will help understanding the 
Very Early Universe properties. 

PACS 11.30.Pb - Supersymmetry. 
PACS 95.35.+d - Dark matter. 

PACS 98.80.Cq - Particle-theory and field-theory models of the early Universe. 



1. — Introduction 

The present and future high energy colliders will hopefully allow the discovery of 
new particles. Many new physics models beyond the Standard Model propose particle 
candidates for dark matter, within the reach of the LHC and ILC. Should a stable, neutral 
and weakly interacting new particle be found, it would be a candidate for dark matter. 
Through particle physics computation of annihilation and co-annihiliation diagrams, it 
is possible to compute the dark matter relic density [1]. This relic density is then often 
compared to the dark matter density deduced from cosmological observations in order to 
constrain new physics parameters. If the computed relic density is compatible with the 
cosmologically inferred dark matter density, the cosmological model will be reinforced by 
this new particle discovery. However, in case of disagreements, two possible paths can 
open: either the new physics particle model is not correctly designed, or the cosmological 
model is missing ingredients, such as quintessence [2] or reheating [3]. In [4, 5], we showed 
that such dark components can modify the computed relic density by several orders of 
magnitudes, and we introduced parametrizations to characterize the Very Early Universe 
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properties. Using this kind of parametrizations, and combining with particle physics data, 
it will therefore be possible to determine the Very Early Universe properties beyond the 
standard cosmological scenario. In the following we will consider Super symmetry (SUSY) 
as an example and discuss the use of relic density to constrain cosmological properties. 

2. — Relic density calculation 

The density number of supersymmetric particles is determined by the Boltzmann 
equation and takes the form: 

(1) ^ = -3Jfn-M(n 2 -n^), 

where n is the number density of supersymmetric particles, (av) is the thermally averaged 
annihilation cross-section, H is the Hubble expansion rate and n eq is the supersymmet- 
ric particle equilibrium number density. The expansion rate H is determined by the 
Friedmann equation: 

(2) H 2 = —(p rad + p D ). 

p ra d is the radiation energy density, which is considered to be dominant before BBN in 
the standard cosmological model. pr> is introduced in Eq. (2) to parametrize the expan- 
sion rate modification [4] and can be interpreted either as an additional energy density 
modifying the expansion (e.g. quintessence), or as an effective energy density which can 
account for other phenomena affecting the expansion rate (e.g. extra-dimensions). 

The entropy evolution can also be altered beyond the standard cosmological model. 
In presence of entropy fluctuations we give the entropy evolution equation: 

(3) ^ = -3Hs + E D , 

at 

where s is the total entropy density. in the above equation parametrizes effective en- 
tropy fluctuations due to unknown properties of the Early Universe, and is temperature- 
dependent. 

The parameters pn and Y>d can be regarded as independent. Entropy and energy 
alterations are considered here as effective effects, and can be generated by curvature, 
phase transitions, extra-dimensions, or other phenomena in the Early Universe. In a 
specific physical scenario, these parameters can be related, as for example in reheating 
models [3]. 

The radiation energy and entropy densities can be written as usual: 

7T 2 27T 2 

(4) Prad = 9eff CO — T 4 , S rad = h eff (T) — T 3 . 

We split the total entropy density into two parts: radiation entropy density and effective 
dark entropy density, s = s rac i + sp. Using Eq. (3) the relation between sp and can 
then be derived: 



(5) E D = J——y/i + p D T 
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with 

(6) gl/> = 2!^( 1+ T*£* 

V<kff V 3 h eff dT 

Following the standard relic density calculation method [1], Y = n/s is introduced, and 
Eq. (1) yields 



dY _ misp I 7T 1/2 ^ 1 + s D 



{™){Y*-Y? q ) + - ™° 



(h eff {T)%T*) (l + s D ) 



dx x 2 V 45G a VVl + Pr- 

(7) 

where x = mi sp /T, m^ sp is the mass of the lightest supersymmetric relic particle, and 



sd ~ __ Pd 



(8) ^ %(T)fT3 ' ^-^TfTi' 

and 

with z running over all supersymmetric particles of mass ra$ and with ^ degrees of 
freedom. Integrating Eq. (7), the relic density can then be calculated using: 

(10) m 2 = mis P s f r ° h2 = 2.755 x 10 8 F m Jsp /GeV , 

Pc 

where the subscript refers to the present value of the parameters. In the limit where 
Pd = <§ j d = H j d = 0, standard relations are retrieved. Using Eqs. (1-10) the relic density 
in presence of a modified expansion rate and of entropy fluctuations can be computed 
provided pr> and sb are specified. For pjj we follow the parametrization introduced in 
Ref. [4]: 

(11) Pd = K P Prad(T BBN )(T/T BBN ) np , 

where T BB n is the BBN temperature. Different values of n p leads to different behaviors 
of the effective density. For example, n p = 4 corresponds to a radiation behavior, n p = 6 
to a quintessence behavior, and n p > 6 to the behavior of a decaying scalar field. k p is 
the ratio of the effective energy density to the radiation energy density at BBN time and 
can be negative. The role of p B is to increase the expansion rate for p B > 0, leading to 
an early decoupling and a higher relic density, or to decrease it for p B < 0, leading to a 
late decoupling and to a smaller relic density. To model the entropy perturbations, we 
follow the parametrization introduced in Ref. [5] : 



(12) s D =n s s rad (T BBN )(T/T 1 



BBN 



This parametrization finds its roots in the first law of thermodynamics, where energy 
and entropy are directly related and therefore the entropy parametrization can be similar 
to the energy parametrization. As for the energy density, different values of n s lead to 
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different behaviors of the entropy density: n s = 3 corresponds to a radiation behavior, 
n s = 4 appears in dark energy models, n s ~ 1 in reheating models, and other values 
can be generated by curvature, scalar fields or extra-dimension effects. k s is the ratio of 
the effective entropy density to the radiation entropy density at BBN time and can be 
negative. The role of sd is to increase the temperature at which the radiation dominates 
for sd > 0, leading to a decreased relic density, or to decrease this temperature for < 0, 
increasing the relic density. Constraints on the cosmological entropy in reheating models 
have been derived in [6]. 

3. - SUSY constraints 

We now consider the effects of the parametrizations described in the previous section 
on the supersymmetric constraints. Using the latest WMAP data [7] with an additional 
10% theoretical uncertainty on the relic density calculation, we give the following favored 
interval at 95% C.L.: 

(13) 0.088 < n DM h 2 < 0.123 . 
The older dark matter interval is also considered: 

(14) 0.1 < Vt DM h 2 < 0.3 . 

One million random SUSY points in the NUHM parameter plane (/i,m^) with mo = 
rai/2 = 1 TeV, A = 0, tan/3 = 40 are generated using SOFTSUSY v2.0.18 [8], and 
for each point we compute flavor physics observables, direct limits and the relic density 
with Superlso Relic v2.7 [9, 10]. In Fig. 1, the effects of the cosmological models on 
the relic density constraints are demonstrated. The first plot is given as a reference for 
the standard cosmological model, showing the tiny strips corresponding to the regions 
favored by the relic density constraint. In the second plot, generated in a Universe with 
an additional energy density with k p — 10 -4 and n p = 6, the relic density favored strips 
are reduced, since the calculated relic densities are decreased in comparison to the relic 
densities computed in the standard scenario. The next plots demonstrate the influence 
of an additional entropy density compatible with BBN constraints. The favored strips 
are this time enlarged and moved towards the outside of the plot. This effect is due to a 
decrease in the relic density. These figures show that a modification in the cosmological 
scenario can completely modify the calculated relic density and lead to different shapes 
of the favored parameter regions. 

4. — Inverse Problem 

The determination of the supersymmetric parameters using particle physics observ- 
ables can on the other hand give access to global properties of the relic particle decou- 
pling period [11]. Let us consider the NUHM example point (mo = = 1 TeV, 
rriA = /i = 500 GeV, Ao = 0, tan/3 = 40), which gives a relic density of Qh 2 ~ 0.11 in 
the cosmological standard model, compatible with WMAP results. The effects due to 
the presence of effective energy or entropy in the Early Universe are presented in Fig. 
2: the first plot shows the influence of an additional effective density on the computed 
relic density. We note that when k p and n p increase, the relic density increases up to a 
factor of 10 5 . The second plot illustrates the effect of an additional entropy density, in 
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Fig. 1. - Constraints on the NUHM parameter plane (/i,raA), in the standard cosmological model 
(top left), in presence of a tiny energy overdensity with k p — 10 -4 and n p — 6 (top right), and 
of an entropy overdensity with k s — 10 -5 and n s — 5 (bottom left), with k s = 10 -4 and n s — 5 
(bottom right). The red points are excluded by the isospin asymmetry of B — >• K*j, the gray 
points by direct collider limits, the yellow zones involve tachyonic particles, and the dark and 
light blue strips are favored by the dark matter constraints of Eqs. (13) and (14) respectively. 



absence of additional energy density. Here when k s and n s increase, the relic density is 
strongly decreased down to a factor of 10 -14 . The dark lines delimit the zones which are 
compatible with WMAP data. The determination of the NUHM parameters provides 
therefore interesting constraints on the cosmological properties of the Early Universe. 

It is important to point out that all the cosmological scenarios previously described 
are equivalent from the point of view of the cosmological observations: there is no way 
to distinguish between them with the current cosmological data. 

5. — Conclusion 

We discussed the use of the LHC and ILC data, as well as flavor physics constraints, 
together with the dark matter relic density evaluation, to explore the properties of the 
Very Early Universe. In particular, we provided par ametrizat ions of the entropy content 
and the expansion rate for this purpose. We have shown that a disagreement between the 
computed relic density and the observed dark matter density may be a sign of deviation 
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Fig. 2. - Influence of the presence of an effective energy density with n p = 6 (left), and an 
effective entropy with n s — 5 (right). The colors correspond to different values of Qh 2 . The 
black lines delimit the regions favored by WMAP. The favored zones are the lower left corners. 



from the standard cosmological model. Should high energy colliders discover candidates 
for dark matter, such studies of the Early Universe properties will become possible. 
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